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CALIBRATION OF TRANSITION MATRICES 


Section 5.3.2 discusses the way that the current 
version of the Activity Block is calibrated. The regression 
equations are simply adjusted so as to yield a dependent variable 
identical with the most recent observations. This method, 
although simple and straightforward, May not always be con- 
venient. This is because the transition probabilities that 
one needs may not exist, or because they may be inaccurate. 
Transition probabilities are far more difficult to obtain 
than state distributions, and are often approximated from 


more than one source, and are hence not wholly reliable. 


The state distribution method discussed below is 
an attempt to adjust the probabilities so as to generate a 
desired distribution exactly. It is assumed that the original 
Matrix is for the most part quite “qood” anc tha at sneule 
be changed as "little" as possible. The problem then is to 
find the smallest changes possible in the a so as to 
generate the measured distributions exactly, and so as to 


ensure that each row of the given matrix still sums to one. 


This no doubt sounds very vague. We can be more 


specific by stating the problem mathematically. 
Mathematical Formulation of the Problem 


(a) Let p;. be an element of the original transition 
13 


matrix and let x, 


a : +O. De. th orrespondi 
a Ps e ec sponding 


j ij 

element in the adjusted transition matrix. We 

wish to find ae SO aie ale is as "small as 
possible" for all i and j and so that the resulting 


X matrix will yield no deviations from the historically 


observed distributions. 
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(b) More rigorously, we are given: 


(i) 4 Pies Pi3 = Original transition matrix. 
(ii) S(t) = a vector equal to the number of persons 
in the three employment states at 


time t. 


That is; S'(t) = (S53 (t) 55 4c) eo te 


where S. (t) = number of people in state i 


at time t: 


aia), | Oo" (O1, G2, G3) 


a vector equal to the actual number of 


i] 


people in the three states at time t+l. 


(iv) S(e) 2 = S(t+1) e456 


That is; the actual distribution vaceci me 
t+l1 does not equal the simulated distribution 


at time t. 


The problem now is to find: 


X = x.. = a transition matrix such that S' (t)X = o 
13 


and such that X is as close as possible to P. 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Toronto 


https://archive.org/details/31/61117013516 


(c) We will define the degree of "closeness" to be 


n n 2 
J» > py (ie, Cel eoe es 
ial hake od eine 


That is, we wish to make the sum of the squares of 
the deviations between the elements of the two 


Matrices as small as possible. 


Ca)" “LE we define Ps and x. to be column i of the 
respective matrices, we can formulate the problem 


precisely as follows: 


n 
ee c 2 : - 
Minimize J = 2 (x, P,) (x, P,) 
i=l 
subject to S'X; = oF ies al eg & 
and n 
> x, =1 
isl 


The first constraint forces the simulated distribu- 
tion to equal the actual distribution, while the second 


constraint forces each row of the derived matrix to sum to: i. 


Solution 


The Lagrangean is, 


n n 

L = 0 <2" ee = ]) = ats (S x, — 0 ) 
and 
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or Ex = l 
1 oL 
Sant | SO Se. = Oro 
ae) " a A 
’ = — 
or S xX, = Os isl, n 
From (1) we have 
Lx; = tps on + LAS 
or le=tl1l+nprt LASS 
or a rA.S 
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Substituting (4) in (1) gives 
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Substituting (6) in (5) gives 


=i = Pe + (9,-S'p.)s (7) 
Ss‘s 
n 
or Soa oo Oy _-O. 5 
ij Pi; bP cr (8) 
n 
LS. 
kel * 


From which, 


Equation (8) can now be solved for ae ina 
straightforward manner, since all the variables on the right 
hand side are given. The equation shows that the change in 
the transition probability between state i and state j is 
proportional to the total deviation between the simulated 
total and the actual total in state j (the bracketed term), 
proportional to the number of people in state i at time t, 
and inversely: proportional to the sum of squares of the 
totals in each state at time t. In writing the computer 
program to calibrate the matrices it will be necessary to 
test for the non-negativeness of Hous Since it can be seen 


that non-negativeness is not guaranteed by equation (8). 


This method, it should be noted, has not been 


tested in practice. It is only offered as a first attempt 


to come up with some reasonable method of calibrating tran- 


sition matrices. There is no way of telling, at the moment, 


how satisfactory this method might turn out to be in practice. 
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APPENDIX D.2 


VALIDATION OF THE ACTIVITY STATUS BLOCK 


VALIDATION OF THE ACTIVITY STATUS BLOCK 
et TALUS BLOCK 


The parameters of the Activity Block were validated 
by simulating the period from April 1972 througniaprit 1o7o. 
The results of these validation simulations are contained in 
the attached graphs and tables*. The discussion below 


pertains mainly to the graphs. The tables are simply a more 


detailed representation of the same phenomena. 


The graphs attempt to show two things. The first 
is the degree to which the simulation adequately describes 
the unemployment behavior of certain groups in the population. 
To show this, the number of simulated unemployed persons in 
a particular region, or a particular age-sex group, or a 
particular marital status are compared with the actual 
number of unemployed in that group. The second thing the 
graphs attempt to demonstrate is the effects of the adjustments 
that were made to the original parameters. It will be 
recalled that three adjustments were made. The first took 
cognizance of the fact that POLSIM dichotimizes the labor 
force into two groups: Class A individuals who never become 
unemployed, and Class B individuals who are subject to 
unemployment. Adjustment was necessary to allow for this 
since the original parameters were derived on the assumption 
that anyone in the entire labor force could become unemployed. 
The second adjustment calibrated the age-sex parameters to 
account for structural changes in the economy that were 
assumed to have obtained after the original parameters were 
estimated. The third adjustment took account of the fact 
that one's labor force behavior depends on his marital 
status and region. Each age-sex matrix in the original 
eee was transformed into 10 marital Ape as 
ericee: each reflecting the adjustment necessary for the 
particular class. 


Ce Ee 
Because of the extensive nature of the validation tables, they 
included in the present report. They are available on 
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The graphs can be divided into 18 sets of four. 
Each set plots the number of unemployed in a particular 
group (total population, 5 regions, 2 marital status. groups, 
8 age-sex groups, and 2 sex groups) for each of the 13 months 
simulated. Within each set the following four comparisons 
are made, each on a separate graph: Original simulation 
versus the final (all 3 adjustments) simulation; the original 
simulation versus the adjustment for Class A persons; the 
Class A adjustment versus the calibrated Class A adjustment; 
and the calibrated Class A adjustment versus the calibrated 
Class A adjustment with disaggregation. As a common reference, 
the actual values for the particular group are also plotted on 
each graph. It is thus possible to examine the cumulative 
and net effect of the three adjustments as they affect any 


given group. There are 76 graphs in all (4x18). 


In examining the graphs (or the tables) it is clear 
that only the final (all 3 adjustments) version adequately 
simulates the population over all relevant variables (age, 
sex, region, and marital status). The other three versions 
all contain significant errors when distributions over marital 
status and region are examined. This is as one would expect, 
since the other three versions are stratified only on age and 
sex. If ne calibrated and final versions are compared for 
a given age-sex group on the other hand, it can be seen that 
there is little to choose between them. Again this is as 


one would expect, since in this case both depend upon age 


and sex. 
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The adjustment for Class A persons performed very 
well, as can be seen by comparing it with the original 
Simulation. The two should be identical, and for all practical 
purposes, they are. It can be seen that the definition of 
Class A individuals is not quite perfect. We tend to allocate 
too many people to this class, as evidenced by the drop in 
total unemployed in the base month. The simulation tends 


to correct for this, however, and in general performs more 


than adequately. 


The size of the errors are well within the tolerance 
range one would expect for this kind of simulation. The 
error in the total unemployed in the final version ranges 
from 5.6% in the base month to .2% in July. The mean error 
is 2.3%. The error in total employed ranges from .4% in the 
base month to 0.0% in July and August. In no simulated month 


does the error rise above .2%. The mean error is .14%. 
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APPENDIX E 


THE MARKET INCOME BLOCK 


APPENDIX E.1 


_ DEFINITION OF THE PARAMETERS 


DEFINITION OF THE PARAMETERS 


The following is a list of all the parameters that 
are used by the Market Income Block. The constant parameters 
are read in from a single tape and stored in the named Common 


Blocks. The variable parameters are read in from cards. 
4 Constant Parameters 
(a) COMMON/ALPHA/PENS (18,2,9) 
PENS (1I,J,K) is the cumulative probability that a 
person whose final earnings class is I and whose sex is J will 


have his pension equal to K - class percent of final annual 


earnings. 


tas, ly < 2.5K$ 10 7-8K$ 
2 2.5-3K$ ea 8-9K$ 
3-3 .5K$ 12 9-10K$ 
3 .5-4K$ tie: 10-11k$ 
4-4 ,5K$ 14 11-12K$ 
A. 5-5SKp 15 12-13K$ 
5-5 .5K$ 16 13-14K$ 
5 .5-6K$ ung 14-15K$ 
6-7K$ Jigs) 15K$+ 
J = 1, Male 2, Female 
K = 5% 35% qd 65% 
15% 45% 8 75% 
25% 55% 9 90% 
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(b) COMMON/BETA/ANWAGE (14,14 ,3,2) 


ANWAGE (I,J,K,L) is the annual employment income 
transition array in real dollars (1970 dollars). For age class 
K and sex L, ANWAGE (I,J,K,L) is the cumulative pEODabiditcy of 
moving from income class I (where income is in 1970 dollars) 


in year t to income class J (1970 dollars) in year t#l. 


f= Input income class (row of matrix) 
1 1-999 8 7K-7999 
2 1K-1999 9 8K-8999 
3 2K-2999 10 9K-9999 
4 3K-3999 Li 10K-14999 
5 4K-4999 2 15K-19999 
6 5K=5999 13 20K-24999 
Hf 6K-6999 14 25K+ 

J = Output income class (column of matrix) 


_ coded as above 


K = Age class 
1 14-35 
2 36-45 
es 46+ 
i= sex 
AE Male 


2 Female 
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(c) COMMON/GAMMA/WEEKWG (10,10,4,2) 


WEEKWG (I,J,K,L) is the weekly wage rate transition 
peray ini realy (1969) dollars. For: age class K and. sex LL)! WEEKWG 
(I,J,K,L) is the cumulative probability of moving from wage class 
Pawiere the wage rate is stated in 1969 dollars) in year t to 


wage rate class J (1969 dollars) in year t+1. 


I = Input wage class (row of matrix) 
u 0 6 101-120 
2 1-40 a 121-160 
3 41-60 8 161-200 
4 61-80 9 201-240 
5 81-100 10 241+ 

J = Output wage class (column of matrix) 


~ coded as above 


K = Age class 
z 14-24 3 36-45 
2 25-35 4 46+ 
Ls sex 
1 Male 


2 Female 
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(d) COMMON/DELTA/INITY (57,2,10,3,10) 
INITY (1I,J,K,L,M) is the mean full employment income 
in 1971 dollars for a person in age class I, sex J, education 
class K, marital status L, and province M. 


I= Age 


Pe smate a fo. eS LS ot .5. 57 2s 10 


J = Sex 
ih Male 
2 Female 
K = Education class 
1 no school 
Pe some elementary 
3 elementary complete 
a some high school 


5 high school complete 


6 some University 
iL CAAT graduate 

8 Bachelor's level 
9 Master's 
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L = Marital Status 


1 Single 


Ds married 
3 other 

M = Province 
bk NFLD 6 Ont. 
2 ID Aoyagi es Man. 
3 INGO 8 Sask. 
4 2B. 9 Alta. 
5 P.Q% 0: B,C; 


(e) COMMON/EPSLON/G (6,2,6) 


G (I,J,K) is the one year exponential growth factor for 
real incomes for people in age class I, sex J, and education class 
Peeeif YREAL (t) is the mean constant dollar income in year © for 


a given class, then the constant dollar income in year t+n will be 


YREAL (ttn) = YREAL (t) * e9” 


fe ege Group 
if 14-19 4 30-34 
2 20-24 5 35-44 


3 25-29 6 45+ 
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J = Sex 


1 male 
2 female 
K = education class 
1b <gr 9 A some post-secondary 
2 or 9-11 5 Univ. graduate 
a Gels ornaL3 6 post graduate 


. 


(£) COMMON/ZETA/DIVTRN (13,13,4), GNTRAN (13,13,4) 


DZERO. Cl3, 47372) 7 G2BRO. (134732) 


(i) DIVTRN (I,J,K) is the dividend cumulative transition 


matrix for people in age Class K. 


I and J are the dividend classes: 


ak 0 8 3K-4K 
2 1-250 9 AK-5K 
3 251-500 . 10 5K-6K 
4 501-750 ita 6K-7K 
5 751-1000 AP 7K-8K 
6 1K-2K t3 8K+ 

7 2K-3K 

Kas age class 
l 14-35 3 50-65 


2 35-50 4 65+ 
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(ii) GNTRAN (I,J,K) is the interest cumulative transition 
matrix for people in age class K. The indices are 


the same as for DIVTRN. 


(iii) DZERO (I,J,K,L) is the cumulative probability of 
moving from initial dividend class L to dividend 
class I for people in age class J and income 


class K. The indices are as follows: 


final dividend class 


1 0 8  3K-4K 
2 1-250 9  AK-5K 
3 251-500 10) | 56K 
4 501-750 1 6K=7K 
5 751-1K 12  7K-8K 
6 1K-2K 13. 8K+ 

a 2K-3K 


age class 


a 14-35 3 30-65 


) 35-50 4 65+ 


income class 


Rd-lit bd 


i= nwt lal aividend class 


2 = 1-250 


(iv) GZERO (I,J,K,L) is the cumulative probability of 
moving from initial interest class L to interest 
class I for persons in age class J and total income 


class K. The indices are the same as the above. 
es Variable Parameters 
(a) COMMON/ETA/KYR,RINFL 
Caan RY R is the year being simulated (1972, etc.) 
(ii) RINFL is the annual rate of inflation that is 
assumed to obtain from 1972 forward. It 
is expressed in percentage terms. For 


example, 6% inflation would be input as 
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‘APPENDIX E.2 


THE ESTIMATION OF TRANSITION MATRICES 


THE ESTIMATION OF TRANSITION MATRICES 
——— EE EEN SAE RL 


de The most direct method of estimating transition matrices 
is to obtain counts of the number of events that correspond to 
each cell of the matrix we are trying to derive. The probabilities 
are then obtained by simply dividing the number of counts ina 
given cell by the row total Sool esti ones that cells. Usualiv, 
the counts are obtained from a sample of the population for which 
the matrix is to apply. A problem that arises in this procedure 
is that it is not obvious how large the sample should be so as to 
ensure that the sample parameters are a close approximation to the 
corresponding population parameters. This appendix addresses itself 


to this question. 


25 Consider row i of the transition matrix and let aa be 
the probability that a person in that row will also be in column ice 
Then ore = aan will be the probability that he is in any other 
columns. To keep the notation simple we will write P and G50 
ae and qi5 in what follows. 

In a sample of -N people (N is the row total) let x be 
the number of people who fall into cell j. Then the distribution 
of x is binomial. Our problem is to determine how large N should 
be so as to ensure that x approximates the true population para- 


meter, or to ensure that x/N approximates P. 


3. More formally, let r be the sample estimate of P. We 


wish to choose N such that 


< < 
ysl 1 eSB) es sl 


where a is the maximum value of the error that can be tolerated. 
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If N is large and if neither P nor Q is too, closeto 
zero, the binomial distribution can be closely approximated by a 


normal distribution with standardized variable given by 


Er eselN2 
NPq 


That is, Z is normally distributed with zero mean and variance one. 


4. Dividing numerator and denominator by N, we have 
Nl aa 
N 
< < 


If we choose a confidence interval of b, what we wish is that 


< é 
ee > GAL ra 
Prob Pq Pq = b 
N N 


where Z is N(0,1) 


feo is .95 and a is .05 then 


= 
pO Pe eZ 
eq 
N 
or 
N = 42Po¢ 
20025 
fow Pq is a maximum when P = q = .5 


Thus, in general, it is necessary that 
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That is, to ensure that with Probability 2.95 Fis wiethan ene 


range P + .05 it is necessary that there be 400 people in the 


row. 


The number of people in the row depends on the size 
of error that we choose and the confidence interval we decide 


upon. Some of the possibilities are given in the table below. 


Maximum Error 


O01 205 mek 
og 16,600 664 166 
Confidence eS) 10,000 400 100 


Interval - 68 27500 100 25 


Number of people in a row to ensure given 


maximum error with given confidence. 


It should be noted that the N specified above isa 
lower bound if the transition matrix has more than 2 columns. To 
say that we can be 95% confident of the error in a given element 
in a given row is not to say that we can have the same confidence 
in the whole row. The larger the row, the greater the number of 


observations required. 
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APPENDIX E.3 


COMPUTER PROGRAM INCOME 
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